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Abstract: We present a new formulation of the staggered fermion on the D-dimensional 
lattice based on the SO {2D) Clifford algebra, which is naturally present in the action. The 
action of the massless staggered fermion is invariant under the discrete rotation and the 
SO{2D) chiral and other discrete transformations. From transformation properties of the 
fermion, we find two local meson operators (one scalar and one pseudoscalar) in addition 
to two standard meson operators. 
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1. Introduction 

The fermions on lattice are bound to suffer from the doubhng problem. In the early 
stage of the development of the lattice gauge theory, it is proposed that massless modes 
due to the doubling may be regarded as internal (flavor) degrees of freedom ||, ||. The 
naively discretized Dirac fermion was formulated as a staggered fermion in Ref. Q. The 
reconstruction procedure was found to give spinors with flavor degrees of freedom for the 
staggered fermion with one component spinor on each site [^. 

Since the spinor and flavor degrees of freedom originate from its geometrical structure, 
the staggered fermion must transform nontrivially under the rotation. To clarify it is one of 
the main motivations of the present paper. In this paper, we consider the one component 
staggered fermion in D dimensions, which has 2^ degrees of freedom due to the doubling. 
The number coincides with the dimension of the spinor representation of SO{2D). Indeed, 
we will show that the SO {2D) Clifford algebra is naturally present in the theory and it 
plays a crucial role to study symmetries of the staggered Dirac operator. We study the 
(discrete) rotational symmetry, the chiral symmetry and other discrete transformations 
such as parity, charge conjugation and time reversal. 

Symmetries of the staggered fermion have been discussed earlier in Ref.Q. Our ap- 
proach differs from theirs on two important points: 1) The S0{2D) structure is fully 
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respected ; 2) We use a different definition for the rotation, which, we beheve, is more 
suitable for the staggered fermion. The details will be discussed in section 3. 

With the knowledge of symmetries, we find two extra meson operators, more than those 
well-known scalar and pseudoscalar operators. These operators have been overlooked in 
simulations and it is very important to take them into account for a reliable study. 

2. The Staggered Fermion and SO{2D) Clifford Algebra 

In this section, we present a new expression of the staggered Dirac operator based on the 
SO {2D) Clifford algebra to be found for the Z)-dimensional lattice. For the purpose, it is 
important to understand the geometrical structure associated with the staggered fermion 
and its relation to the SO {2D) Clifford algebra. This will be explained in the next subsec- 
tion. Then we proceed to the expression of the staggered Dirac operator that respects the 
algebra. 

2.1 Geometrical structure of the staggered fermion 

To clearly state the geometrical structure, we classify links, plaquettes and hypercubes.* 

Consider the link directed to the positive fj, direction from the site n^. We denote it by 
(n, fi) and assign the sign factor, (—1)"^ . As for the link directed to the negative direction, 
we write it as (n, — /x) and the associated sign factor must be (—1)"''"'"^. The plaquette 
(n, jj.i') is defined as the ordered links. 



On each link of {n^^v)^ we have a sign factor. We refer to the plaquette with (+ -|- +-I-) 

or ( ) sign factors as a cell-plaquette, while that with a set of mixed signs as a 

pipe-plaquette. On any two dimensional surface of the D-dimensional lattice, we find the 
checkered pattern, or the Ichimatsu pattern, formed by the cell- and pipe-plaquettes. 

Note that there are distinctive hypercubes on the lattice: Those formed solely by cell- 
plaquettes with all plus (minus) sign factors will be called microcells (macrocells) by the 
reason to be explained shortly. The coordinates (n^) of sites on a microcell can be written 
as = 2N^ + r^, where A'^^ is some integer and takes the value of or 1. As for a 
macrocell, those are written similarly as = + \ + r^. 

After having introduced the notions of micro- and macrocells, we observe a slightly 
more detailed structure on a two dimensional surface than the Ichimatsu pattern. The 
D = 2 example is shown in Fig.l, where a microcell (macrocell) is shown as a shaded 
(gray) square. We would find the same pattern on any two dimensional surface of the 
D-dimensional lattice. It is noteworthy to realize that the pattern is invariant under 
the translation by twice the lattice spacing, 2a. We call this translation the modulo 2 
translation. 

*Some notions presented here are refined versions of those reported ear her ffil l7|. 



(n, liv) = (n, /i)(n + //, v){n + fi + v, -fJ-){n + v, —v) {fi < u). 
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After the reconstruction, the Grass- 
mann variables on the sites of a microcell 
are to form a spinor labelled by (A'^)^. 
The relative coordinates (r^) on the Grass- 
mann variables are transformed into spinor 
and flavor indices. In this sense, a mi- 
crocell may be regarded as an internal 
space. The global structm'e such as the 
fermion kinetic term are formed out of 
the remaining geometrical structure. 

The relation of the staggered fermion 
to the S0{2D) Clifford algebra comes from 
the fact that sites on a microcell can be 
regarded as the weight lattice for the spinor 
representation of SO{2D). Therefore the Grassmann variables on sites in a microcell are in 
the spinor representation. Accordingly, we will find that the Dirac operator has a natural 
expression in terms of gamma matrices associated with this S0{2D). 

In the next section, we will consider the symmetry of the action under the 7r/2 rotation 
around the center of a microcell. One could have considered the 7r/2 rotation around a site. 
Actually, the latter rotation was studied earlier in Ref. ||5| and the rotational symmetry of 
the staggered fermion was established. These two rotations are to end up with the same 
rotational symmetry in the continuum limit. We choose the former definition to respect 
the geometrical structure explained above that is inherent to the staggered fermion. ^ 

2.2 The staggered Dirac operator 

The single component staggered Dirac operator is given as 

[Dst)n,n' = 2^Vfi{n) — (2.2) 

A* 

where rifj,{n) = r/^(r) = {—l)^''<i^^''. The link variable acts on the fermion variable in 
the representation "rep" of the gauge group. For simplicity, we drop the superscript "rep" 
from the link variable in the rest of the paper. 

In preparation to rewrite the operator in eq. ( |2.2| ) to act on the coordinates (r, A^) , 
we introduce notations for gamma matrices and their products for the S0{2D) Clifford 
algebra. By e, we denote a bit-valued Z)-dimensional vector: That is, each entry takes the 
value of or 1. For each e, we define the matrix 



^In this sense, (N^) can be regarded as the coordinate of a microcell. 

■'■The it/ 2 rotations around a center of any hypercube other than a microcell also leave the geometrical 
structure invariant. However, we do not further consider these possibilities for the rotational symmetry in 
this paper. 
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Figure 1: The microcells and macrocells 
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(r/.,e-)(n,ri)...{r^,ri,) = ((^^3^^ ® ' ' ' C3 <J^'°) (2.3) 

^ \ 
X (0-3 (g) ■ ■ ■ ® f73® CTi (g)l ® • • • (g) 1) 

/ (ri,ri)---(rD,r^) 

As special cases, we obtain a basis for the SO{2D) Clifford algebra, 

> = r^,0' > = -^r^,e^> (2-4) 



with being the unit vector along the /^-direction. Note that 7^ in eq. ( |2.4D are for an 
irreducible representation of the SO {2D) Clifford algebra while they are reducible as for 
the SO{D) Clifford algebra. In terms of the basis given in eq. (^7^), eq. ( p. 3D may be 
expressed as 

(r^,e)(ri,r!)-(rB,r' ) = ((niTl)^' •■■ (nD7D)^''7M), , (2-5) 

We may rewrite (-Dst)n,n' as an operator acting on the coordinates {r,N). That is 
realized by using the expression (^^), 



{Dst)n,n' = {D St) {r,N),(:r',N') — "^^(^ ^i,e) {r,r'){D p,^e){N ,N')- (2-6) 



Here D^^g is the generalized difference operator. 



{D^,e){N,N') = TTn y^X~^y'''i^V)(N,N')i (2-7) 



2D 

s 



with 

< 



(Vp(Ar,Ar') 



(2.8) 



Note that cr^ = and 1 give the backward and forward difference operators respectively. 
As observed in eq. ( |2.7D , the D-dimensional bit- valued vector a is dual to e. 

In the last expression of eq. ( |2.6| ), the Dirac operator is decomposed into the difference 
operator acting on the microcell coordinate and the matrix acting on the relative 

coordinate (r^). Accordingly, we collectively treat the fermion variables on sites in a 
microcell by introducing the new variables: 

^'r(iV) = 6iV+r, ^r(iV) = ^TV+r- (2.9) 

The new fermion variables ^'r(-^) and ^r(A^) belong to spinor representations of SO {2D). 
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3. Symmetries of Staggered Fermion 



The single component staggered fermion action is invariant under the discrete rotation, the 
chiral transformation, parity and charge conjugation. Here we explain these symmetries. 
In particular, we will find that the S0{2D) Clifford algebra plays a vital role to describe the 
rotational discrete symmetry. Also, the well-known chiral transformation will be identified 
with the U{1)a chiral transformation associated with the SO{2D) Clifford algebra. 



3.1 Rotation 

Let us see how the coordinates are transformed under the 7r/2 rotation around the center 
of the microcell. Without losing the generality, we may take the microcell attached to 
the origin and rotate the system around the center of the microcell. We consider the 
/x(<)z/ rotation. It is easy to find out the following transformation rules: For the relative 
coordinates. 



' p, 



(3.1) 



p = v- 



and for the coordinate Np, 



N. 



p, P^P,^, 
-Nu, p = p, 



(3.2) 



p: 



p = v. 



Here the subscript R implies rotated quantities. For the simplicity of the notation, we do 
not write the subscripts pv explicitly. 

The transformation of the generalized difference operator is found to be 



i-iy D,^,^^{Nr, p = p, 
{-ir^D,,,^iNR,N'j,), p = v, 



(3.3) 



where 



£p, p^ p,iy, 

£u, P = P, 



(3.4) 



P 



U. 



The staggered Dirac operator consists of two factors, Dp g and Tp^g. Having obtained 
the transformation property of the former, we need to find the transformation of the latter 
that compensates the sign factors in eq. (|3.3|) so that the Dirac operator transforms as 



= Vp,D,,V^p, 



(3.5) 
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with a matrix V^^- When this is reahzed, the transformations of ^ and ^ are to be 



in order for the single staggered fermion action, 



to be invariant. 

In order to reahze ( |3.5| ) , we require that Fp^^ transform as 



(-ir-r^,,-.^, p = v. 



(3.6) 



(3.7) 



(3.8) 



The matrix V^^ is to found from this condition ( |3.q ). After some efforts, we find the matrix 
of desired property, 



V, 



1 



(3.9) 



Note that the condition (|3.8| ) determines the matrix V^y up to a phase factor. 

By combining eqs. (|3.3|) and (|3.8D, we easily see that the Dirac operator transforms 



as eq. (^^) . Therefore, we conclude that the single staggered fermion action is invariant 
under the discrete rotational transformation. 

it is easy to find the action of V^^ on the basis of the SO {2D) Clifford 



From eq. ( |3. 
algebra: 



7p, P^fJ-,!^, 

l^ii P = ^i 



^ Ip, P^ P-, V: 



lu, P 

t 7m) P 



P, 

V. 



(3.10) 



Here, note that 7^ do not transform as a vector: The action of V^i^ simply exchanges 7^ 
for >. 

Related to the transformation in eq. (|3l0| ), there is a further curious feature of the 
matrix V^p. We studied the -7r/2 rotation around the center of a microcell which leaves the 
system invariant. This intuitively obvious rotation causes a quite nontrivial transformation 
among sites in a microcell, or the weight lattice of the S0{2D) spinor representation. We 
found that this was achieved by the matrix Vfj,^- In eq. ( |3.1C| ), note that negative is the 
determinant of the transformation matrix realized on the matrices 7^ and 7^. Therefore 
the transformation generated on the 2D gamma matrices belongs to 0{2D) but not to 
SO{2D). In the continuum limit, the staggered Dirac operator becomes flavor blind. As 



a result, the transformation in eq. (3^) reduces to the 7r/2 rotation in the Z?-dimensional 
space. 
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3.2 Chiral symmetry 

It is known that the staggered fermion action is invariant under the transformation 

C = (-l)l"ICn, = (-l)l"l+'Cn. (3.11) 

This symmetry forbids the fermion mass term, even for odd dimensions. In this sense, 
the symmetry has been regarded as a sort of chiral symmetry for the staggered fermion, 
though the origin of the symmetry has not been clearly understood. Here we show that 
this symmetry is nothing but the U{1)a chiral transformation associated with the S0{2D) 
Clifford algebra, 

= exp(i0r2D+i)^, "f' = ^'exp(i6'r2D+i), (3.12) 

with 

n -D(-D-i) 

'^2D+i{r,r) = {-i) (-1) 2 (71 • • -7071 • • •7D)(?^,'r ) 

= (cJ3(8)---®fT3)(r,r') = {-lp6r,r'. (3.13) 



Taking 9 = 7r/2 in eq.( |3T2D , we have ^'^(iV) = f(-l)l''l^^(7V) and ^^(A^) = i(-l)l''l^'r(iV), 
that can be rewritten as eq. ( 3.11| ) for ^ and ^ variables. 

By using the properties, 

{r2D+l,7/^} = {r2D+l,7M} = {r2D+l,r^,e} = 0, (3.14) 

we obtain the relation 

exp (i6'r2D+i)-C'st exp {iOT2D+i) = D^t- (3.15) 



Therefore, the action ( |3.7D is invariant under the chiral transformation in any even as well 
as odd dimensions. 

In the rest of this section, we describe discrete symmetries starting with the parity 
symmetry. 

3.3 Parity and charge conjugation 

We take the reflection as the definition of the parity. The reflection along /u direction, P^, 
only affects the /U-components of coordinates, 

and a link variable transforms as, 

[U2N+cf,p\P, 
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Therefore, we find 

The matrix VJ^ = T2d+i^ij. produces the appropriate sign factors by acting on Tp^e- 

From eqs. ( |3.16| ) and ( |3.17|) , we find the Dirac operator transforms as [Dst]p^ = Vp^DstVjl. 
This imphes the action invariance under the parity. In particular, under the simultaneous 
reflections of all the directions, coordinates change their signs, [n^JpT = —n-fi, and the 
Dirac operator is transformed as [Dst]pT = 7i • • • 7D^st7D ■ • • 7i- 

The charge conjugation invariance of the action can be expressed as the following 
condition on the Dirac operator, 

C-^Dst{N, N')C = D,t{N', N). (3.18) 

First we note that the generalized difference operator, D^^^{N, N'), produces a sign 
factor under the exchange of and A^': 

Dp,e{N',N) = {-1Y^+^D^4N,N'). 



This is due to the property of VJ^ given in eq. (2^) under the exchange of and A^', 

\\^fi)iN',N) - (-1) (V^ )iN,N')- 



To obtain eq. ( 3.18| ), F^^^ has to produce the same sign factor under the action of the 
charge conjugation matrix 

This determines how the charge conjugation matrix C acts on 7^ and 7^: 

C-\C = 7?'7* = r?'7M, 

C~%C = -r^% = r^'% (3.19) 

with r]' = —1. Therefore our charge conjugation matrix is C = F2d+i-^ 

By combing all the results stated above, we now know that the action of the staggered 
fermion is invariant under the lattice version of the CPT transformation. 

As we have already mentioned, the symmetries of the staggered fermion have been 
studied earlier in Ref. The authors refer to the symmetries as shift symmetry (the 
translation by twice of the lattice constant), rotational symmetry, axis reversal (parity) 



'^The matrix C = 1 satisfies the condition (p.l9|) with rj' — 1 
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symmetry, interchange (charge conjugation) symmetry, U{l)v (fermion number) symmetry 
and U{1)a (chiral) symmetry. As the spatial rotation, they chose a site as the center of 
the rotation, while we did the center of a microcell (hypercube) 0. We have chosen the 
center to respect the geometrical structure of the staggered fermion. It may be interesting 
to realize that two rotations are related by a dual transformation: The center of a microcell 
is a site on the dual lattice. Another major difference of our approach from that of Ref. Q 
is that we utilized the S0{2D) Clifford algebra to express all the symmetries. 



4. Scalar Operators 

The results we have obtained help us to write down operators with definite rotational and 
discrete symmetry properties. Let us consider the scalar mesons that can be constructed 
with "local variables" , in the sense that all the fields are associated with a single microcell. 
We find four such operators as given bellow: 

r 

M2{N) = Y,^r{N){TD+l)r,r'^r'{N), 

MsiN) = J2 ^r{N){%)r,r'^r'{N), (4.1) 



where 



for D = even and 



r^+i = (72 di • • • £72 (8) CTi = {-i) 2 7i • • • 7i? 



Td+1 = (7i (g) £72 ■ • • O O"! = ("^ ^ 71 • • • 7^ 

for D = odd. The presence of r/j+i for odd dimensions is due to the fact that the spinor 
representation is reducible as a representation of the SO{D) algebra. The operators M3 
and M4 are scalars since 7^ are simply exchanged under rotations, as we have observed 
earlier in eq. ( 3.10| ). 



It is quite important to realize that only Mi and M2 are considered in studies of the 
staggered fermion. Though the operators Mi and M3 (or M2 and M4) are distinguished 
by the chiral symmetry, they mix up due to the finite mass term. So we have to resolve 
this mixing before taking the chiral limit. 



5. Discussions 



We studied exact symmetries of the staggered fermion that are present even before taking 
the continuum limit. The sites on the D-dimensional hypercube, the microcell, form the 
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weight lattice of the spinor representation of SO {2D). That is the very reason why we 
obtain the spinor representation of the rotational group SO{D). We described properties 
of the Dirac operator under the discrete rotation and showed the invariance of the action. 
We also defined the lattice versions of P, PT and C transformations and showed the 
action invariance. It is worth pointing out that field transformations under the rotation 
and discrete transformations are uniquely defined since the fermion variables are in an 
irreducible representation of the SO {2D) Clifford algebra. From our formulation based on 
the SO{2D) Clifford algebra, we obtained the following results. 1) The transformation in 
eq. ( |3.1lD , that had been regarded as a sort of chiral symmetry, is now identified as the chiral 
symmetry associated with the SO {2D) Clifford algebra. 2) Based on properties of fields 
under exact symmetries, we constructed meson operators, including two new operators, 
one scalar and one pseudoscalar operators. This fact implies operator mixings. So we have 
to take account of mixings to reach a reliable result in simulations. 

The staggered fermion, by its construction, has its spinor and fiavor degrees of freedom 
in the geometrical structure. In the ^ variables, the action is invariant under the site-wise 
gauge transformation. After rewriting it in terms of the \I' variables defined in eq. (|2.90 , we 
still have the same gauge symmetry. However, a part of the gauge symmetry distinguishes 
the components of ^, the spinor and/or fiavor indices. This looks peculiar and we expect 
this part of the gauge symmetry is absent in the proper continuum limit. The question 
how it happens certainly deserves a further study. 

The results reported here could be useful for realizing supersymmetry on lattice. In 
our works [^, 0] aiming at a realization of lattice supersymmetry, we presented models 
that possess a Grassmannian symmetry relating bosonic and fermionic variables. The 
geometrical structure of the staggered fermion, explained in the present paper, played a 
crucial role to realize the symmetry. Since we obtained a susy-like transformation for 
the fermion in the naive continuum limit, we expect that the Grassmannian symmetry is 
related to some supersymmetry in an appropriate continuum limit. 

As we used the staggered fermion in (^,0! the expected supersymmetry may be an 
extended supersymmetry. If it is so, scalar degrees of freedom must be hidden in the 
bosonic variables. Therefore it is very important to fully understand rotational properties 
of the theory including the link variables. 
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